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Many astrophysical phenomena exhibit relativistic radiative flows. While velocities in excess of 
v ~ 0.1c can occur in these systems, it has been common practice to approximate radiative transfer 
to 0(v/c). In the case of neutrino transport in core-collapse supernovae, this approximation gives 
rise to an inconsistency between the lepton number transfer and lab frame energy transfer, which 
have different 0(v/c) limits. A solution used in spherically symmetric 0(v/c) simulations has 
been to retain, for energy accounting purposes, the 0(v 2 /c 2 ) terms in the lab frame energy transfer 
equation that arise from the 0(v/c) neutrino number transport equation. Avoiding the proliferation 
of such "extra" 0(v 2 /c 2 ) terms in the absence of spherical symmetry motivates a special relativistic 
formalism, which we exhibit in coordinates sufficiently general to encompass Cartesian, spherical, 
and cylindrical coordinate systems. 
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I. MOTIVATION 

The requirement of energy- and angle-dependent neu- 
trino radiative transfer makes a core-collapse supernova 
a seven-dimensional problem: In addition to the four di- 
mensions of spacetime, the neutrino distribution func- 
tions depend on a three-dimensional momentum space. 
Energy- and angle-dependent neutrino transport has 
been employed in spherically symmetric simulations by a 
number of groups |T|. HJ. [J] In the spatially multidimen- 
sional case, energy- and angle-dependent transport has 
been implemented in an approximate 'ray-by-ray' fash- 
ion that, except for advection, ignores lateral transport 
|3J . Efforts aimed at unapproximated versions are under- 
way as well 0, . For more information on the need for 
detailed neutrino transport, and the history of approxi- 
mate treatments and resulting insights into the explosion 
mechanism, see for example the reviews in Refs. Q>|g. 

With one exception in spherical symmetry all of 
this work has used equations of neutrino radiative trans- 
fer approximated to 0{v/c); but this gives rise to an 
inconsistency associated with energy conservation, a cru- 
cial check of a simulation's accuracy that has implications 
for claims about the explosion mechanism. The diffi- 
culty is that energy conservation can only be assessed in 
the "lab frame" (that is, by an inertial observer, neces- 
sarily at rest at infinity in the general relativistic case) 
[ToL [Tlj . but the 0(v/c) limit of the lab- frame energy 
transfer equation is not equivalent to the 0(v/c) limit of 
the energy transfer equation in a reference frame comov- 
ing with the background fluid 0. To see this, note first 
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that lepton number conservation is based on the Boltz- 
mann equation, which expresses the transport of neutrino 
number: 

L[/] = q/], (i) 

where L[/] is the Liouvill e op erator and C[/] is the col- 
lision integral [13 . fl3l ITiL Il5| . A commonly used energy 
transfer equation is Eq. Q multiplied by e, the neutrino 
energy measured in the comoving frame: 

eL[/] = eC[/]. (2) 

However, fully conservative forms are obtained in flat 
spacetime [l5| (and even in general relativity, in the case 
of spherical symmetry [ToL HHp by multiplying Eq. (JTJ 
by the particle energy measured by a lab frame observer. 
To 0(v/c), this lab frame energy is e(l + hiVi), where 
e hi are the components of particle momentum in the co- 
moving frame, and v\ are the fluid velocity components 
measured in the lab frame. Hence the lab frame energy 
transfer equation is 

e(l + Mi)L[/] = e(l+n i S i )C[/]. (3) 

When / is taken to depend on neutrino energy and angles 
measured in the comoving frame, L[/] contains velocity- 
dependent terms representing Doppler shifts and angu- 
lar aberrations. But the 0(v/c) limit of L[/] gives rise 
to terms of 0(v 2 /c 2 ) on the left-hand side of Eq. J3J|; 
dropping these yields an 0(v/c) lab frame energy trans- 
fer equation that is not equivalent to the 0(v/c) limit of 
Eqs. P and 0. 

There have been two basic approaches to the prob- 
lem of computing both energy and lepton number trans- 
fer with energy- and angle-dependent neutrino transport, 
but only one group has taken special care with global en- 
ergy conservation. Two groups [r^IItI] have employed a 
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variable Eddington factor approach (which might also be 
called an "iterated moment method" [(J) to solve Eq. (J2J 
for the neutrino energy distributions in the comoving 
frame. The integral of Eq. J2J over momentum space pro- 
vides a heat source to the fluid internal energy, and the 
lcpton number transfer follows easily from the momen- 
tum space integral of Eq. (Q. However, consistency with 
Eq. J2J — which bears on global energy conservation — 
has not been addressed by these groups. In the other 
approach 0, El E E3 (which might be called a "di- 
rect method" @), a complete discretization of Eq. 
is solved directly. Its integral over momentum space 
yields the lepton number transferred to the fluid, and 
the source for fluid internal energy follows from the mo- 
mentum space integral of Eq. (J2J). However, in this case 
special care has been taken to ensure that the discretiza- 
tion of Eq. Q is consistent with Eq. J2J, so that global 
energy conservation is maintained at a level capable of 
resolving with confidence the phenomenon of greatest in- 
terest: an explosion with kinetic energy ~ 10 51 erg. 

This explosion energy of ~ 10 51 erg must be discrim- 
inated against a background of ~ 10 53 erg, the natural 
energy scale of the system. This larger reservoir of grav- 
itational potential energy is released during collapse and 
transmuted into neutrinos, preceding the fractional en- 
ergy return to the fluid in the form of neutrino heating 
that leads eventually to kinetic expansion. This trans- 
mutation is a genuine radiation hydrodynamics problem, 
in which the neutrinos and matter are tightly coupled; 
hence it is difficult to justify an argument that errors of 
~ 10 erg in the global neutrino energy are irrelevant 
to the magnitude of errors in the fluid kinetic energy. 
Maintenance of excursions in global energy to an order 
of magnitude below ~ 10 51 erg corresponds to a precision 
of one part in 10 3 . A typical simulation will be carried 
out over ~ 10 5-6 time steps, which requires that energy 
be conserved systematically to better than one part in 
10 8 ~ 9 per time step. This is a severe requirement, one 
that is very difficult to satisfy in a realistic supernova 
model. 

Such precision is a necessary (but not sufficient) condi- 
tion for confidence that simulation outcomes represent re- 
ality. Even assuming all the relevant physics has been in- 
cluded, it is true that conservation of lepton number and 
energy are no guarantee that a model is correct: It is pos- 
sible that a model could conserve total energy, while not 
accurately partitioning the total energy among kinetic, 
internal, gravitational, and other components. Never- 
theless, a model that does not conserve lepton number 
and energy fails an obvious measure of quality control. 
How should we interpret the prediction of a ~ 10 51 erg 
explosion in a model where the total energy varies during 
the course of the simulation by ~ 10 51 erg or more? 

It is true that in spherical symmetry all three groups 
agree qualitatively on a failure to explode Q, 0, , and 
that two of the groups have been shown to compare 
favorably in quantitative terms during some simulation 
epochs 20J; but the need for a careful treatment of neu- 



trino transport that conserves global energy remains nec- 
essary as the field moves towards two and three space 
dimensions. For one thing, assessing the mutual valid- 
ity of the simulations required a comparison with the 
transport treatment that underwent continued improve- 
ment until the monitored global energy was conserved at 
the required level [TT| . Moreover, the simulation com- 
parisons showed some nontrivial differences during the 
outward trajectory and initial stall of the shock. How 
these differences might play out in the multidimensional 
case — which is closer to exploding — remains to be seen. 
The history of qualitative disagreement in the outcome of 
simulations by different groups in multiple space dimen- 
sions p. l2lL 122 . |23| — not to mention the fact that experi- 
ence has demonstrated that small variations in neutrino 
transport treatments can sometimes lead to qualitative 
changes in outcome [24| — indicates that care and caution 
are called for (see also the review in Ref. 8]). 

Because of the distinct O(v/o) limits in energy and 
lepton number transfer, a treatment valid to all orders in 
v/c is arguably simpler. Infall velocities of 0.2c — 0.4c are 
sufficiently high during collapse that the 0(v 2 /c 2 ) terms 
in Eq. J3J are not negligible, and must be retained for 
the purposes of global energy accounting This in- 

dicates, of course, that the 0(v/c) limit is quantitatively 
inadequate, but there is also a practical reason to go fur- 
ther: In the relativistic case, Eq. actually simplifies, 
with the "extra" 0(v 2 /c 2 ) terms disappearing (see Refs. 
[Tol[TlT | for the spherically symmetric case, and Ref. ^3 
and the present work for the spatially multidimensional 
case). In two and three space dimensions, the number of 
"extra" 0(v 2 /c 2 ) terms is greatly multiplied, rendering 
this practical dividend of a relativistic formalism even 
more attractive. Due to the absence of gravitational ra- 
diation in spherical symmetry, general relativistic simula- 
tions pill l| are not qualitatively more difficult than those 
done in the Newtonian+C(u/c) limit. While supernova 
simulations are not yet of sufficient maturity to take on 
the challenge of full numerical general relativity in two 
and three space dimensions — and will thus have a quanti- 
tative deficiency at some level — a special relativistic for- 
malism nevertheless avoids the "extra" 0(v 2 /c 2 ) terms. 
Enhancements of Newtonian gravity can then capture 

some of the general relativistic effects [nun 12H- 

This paper is organized as follows. In Sec. [HJ 
we present the general relativistic Boltzmann transport 
equation and its conservative variants, and describe our 
space and momentum space coordinate choices. Expres- 
sions in the transport equations are specialized to flat 
spacetime in Sec. IIIII in coordinates sufficiently general 
to encompass Cartesian, spherical, and cylindrical coor- 
dinate systems. Section llVI summarizes our results, and 
mentions a couple of other astrophysical examples re- 
quiring relativistic radiation transport. An appendix de- 
tails the relationship between the number- and lab-frame 
energy-conservative equations — a relationship that ought 
to inform the construction of finite-difference representa- 
tions. 
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II. PREPARATION 



We begin by displaying three compact forms of the 
transport equation — valid even in general relativity — 
that involve global spacetime coordinates and momen- 
tum variables reckoned in a frame comoving with the 
fluid with which the transported particles interact. The 
Boltzmann equation for electrically neutral particles is 



"S-rW/gS = C[/l. (4) 

Equivalent in content are the conservative expressions 
for particle number and energy-momentum conservation 
M: 
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which consist of divergences in spacetime and momentum 
space. 

The quantity at the heart of Eqs. (0J-@) is the scalar 
distribution /, a function of spacetime coordinates x^ 
and three- momentum coordinates u l . (Greek and latin 
letters are spacetime and space indices respectively. Un- 
adorned indices indicate components in a global coordi- 
nate basis, barred indices (~) denote an orthonormal "lab 
frame" basis, hatted indices (") indicate an orthonormal 
basis comoving with the fluid with which the particles 
interact ("comoving frame"), and indices with a tilde (~) 
denote a comoving frame momentum space coordinate 
basis.) The scalar distribution function / gives the num- 
ber of particles dN in an invariant spacetime 3- volume 
element dV and invariant momentum space volume ele- 
ment dP [H: 



(7) 



The quantities x, w, and p are 4- vectors, and p is the spa- 
tial 3- vector portion of p. The unit 4- vector w is timclike, 
and defines the orientation of dV: 



(8) 



where g is the determinant of the metric tensor (taken 
to have signature + + H — ) and ep„ pCT is the Levi-Civita 
alternating symbol (eoi23 = +1). The momentum space 
volume element is 



(9) 



where specialization to momentum variables u l in the 
comoving frame has been made in the second line, and 



(10) 



dp 

du 



du 1 
dp* 



(G) 
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TABLE I: Coordinates, metric functions, and derivatives of 


metric functions for common coordinate systems. 




System x 1 x 2 x 3 a b c 


1 da I db 

a Bx 1 b Bx 1 


1 Be 

ar. Bx 2 


Cartesian x y z 1 1 1 








Spherical r 8 4> r r sin# 


1 1 

r r 


cote 
r 


Cylindrical r z 4> 1 r 1 




r 






is the particle energy measured in the comoving frame 
(m is the particle mass, taken to be zero for the purposes 
of photon or classical neutrino transport). The definition 
of / in Eq. Q makes a convenient connection to nonrel- 
ativistic definitions of the distribution function. For an 
equivalent but more geometric approach, see Ref. [l3T |. 

We now describe in more detail the spacetime and 
momentum space coordinate systems we will use in the 
special relativistic elaboration of Eqs. Q-©- In space- 
time, we have a global coordinate basis. In momen- 
tum space, we invoke orthonormal "lab frame" and "co- 
moving" spacetime bases to obtain the Cartesian spatial 
momentum components p % measured in a comoving or- 
thonormal frame, and then transform these to momen- 
tum space spherical coordinates u l (particle energy, and 
two angles specifying the momentum direction). 

We assume flat spacetime, but in order to accomodate 
curvilinear coordinate systems we label our spacetime 
global coordinates (a; M ) = (x 1 , x 2 , x 3 , f) . (In matrix 
expressions of tensor components, rows and columns are 
ordered 1, 2, 3, 0.) The line element is 



ds 2 = g^dx^dx" , 
where the metric is diagonal: 

(^) = diag[l,a 2 ( a; 1 ),6 2 ( a ; 1 ) C 2 ( a ; 2 ) r 



(11) 



(12) 
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The specialization of the coordinates 3 and met- 

ric functions a, b, c in various coordinate systems are 
given in Table [i] Our spacetime coordinate systems are 
"lab frames," so that the equations we derive are Eule- 
rian. 

The orthonormal "lab frame" is accessed with the 
"tetrad" e^p, which locally transform the metric into the 
Lorentz form t]p,v = diag[l, 1, 1, —1]: 

^ ft & v 9fxf — TJp,t>- (13) 

A simple choice is 

(e%) = diag[l, a-^x 1 ), 6" V) c~ V), !]• (14) 

Because this is diagonal, the inverse simply has the 
diagonal entries inverted. 

The orthonormal frame comoving with the fluid — 
indicated by indices adorned with a hat (") — is related 
to the orthonormal lab frame by a Lorentz boost, which 
satisfies 



where 



A l - 



7 + 1 



ViVj, 



and 



A*o = -yvi, 
A°6 = 7, 



(vif - (v 2 f 



-1/2 



(15) 

(16) 

(17) 
(18) 
(19) 

(20) 



The inverse A^p — the boost from the lab frame to the 
comoving frame — is obtained by taking Vi — > —Vi. The 
bars on the velocities are a reminder that these are the 

3- velocity components in the orthonormal lab frame co- 
ordinate system (they are not spatial components of a 

4- vector) . 

We also define a composite transformation (technically, 
also a tetrad) from the orthonormal comoving frame to 
the global coordinate frame: 



where 



r» am. 



(21) 



(22) 



in momentum space, which for massless particles may 
be expressed 



p = e cos"d = efti, 
p 2 = e sim? cosip e £tj 2 
p 3 = e sin?? sin(/3 = e ti 3 



(23) 
(24) 
(25) 



where the hat on hi indicates that this 3-vector pertains 
to the orthonomal comoving frame. (Note that hi is not 
spatial portion of a 4-vector. With both hi and Vi we 
employ a repeated index summation convention that does 
not require one index up and one index down for these 
non-spacetime-tensor quantities; their indices are always 
down.) 

In summary, the neutrino radiation field is a func- 
tion of the variables (x M ) = far, x 2 , x 3 , i) and (u l ) = 



III. ELABORATION 

With these specific coordinate choices, we are ready to 
specialize the expressions appearing in Eqs. 1] IKil) . Three 
expressions are relatively simple, and one is complicated. 

The simplest are the two determinant expressions in 
Eqs. JHJ and ©. The determinant of the metric tensor 
is g, so from Eq. (fl"2*)) . 



abc. 



(26) 



The Jacobian determinant of the momentum space coor- 
dinate change is 



.1,1 I £ 



e 2 sin i?, 



(27) 



which follows from Eqs. (|23I25|) . 

Also simple are the global coordinate basis momen- 
tum components p^ = D 1 ^p^ appearing in the transport 
equations (the composite transformation C^a is given 
by Eq. From the boost, Eqs. I|16I19|I . and co- 

moving frame momentum components in Eqs. (|l(jp. 
I25| l . the orthonormal lab frame momentum components 
pp- = A^jp' 1 are 



(28) 
(29) 



P 



7 6 (1 + hiVi) . 



7 + 1 



The coordinate basis momentum components p M = e^^p^ 
are 



The inverse D 1 ^ = A^e^ is the transformation from 
the comoving frame to the orthonormal lab frame. 

Finally, the neutrino 4-momentum is described 
in terms of its comoving frame components, 

(p^ = {P^iP^iP^iP^j ■ Only three momentum 
variables are independent; we choose polar coordinates 



P 1 


= 


= P\ 


P 2 


= e 2 ^ 


1 2 

= -p , 

a 


P 3 


= e 3 ^ 


1 3 

= v/ 


p° 


= e -^ 


= p\ 



(30) 
(31) 
(32) 
(33) 
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where we have used Eq. (|14|) . 

The final expression to specialize from Eqs. (1 IKil) is 
(du 1 /dp 3 )T 3 fiop^p" , which, being complicated, consti- 
tutes the computational burden of this paper. The first 
factor can be obtained by computing the Jacobian ma- 
trix (dp 1 /du 3 ) from Eqs. (|23fl - l|25|) and taking the matrix 
inverse; the result is 



( du 1 
[dpi 



I cosi? sini^cos^ sin sin <p \ 
sin d cos $ cos <p cos $ sin ip 



V 



o 



e 

sin^ 



e 

cos (fi 
e sin?9 



(34) 



The connection coefficients in the orthonormal comoving 
basis are 



dC^~ 

i- n>- v>- pi- vp t i- y,>- P g ■ 



(35) 



where 



1 ll<7 ( d 9av . dgap dg Vf . 



vp 



\ dx p dx v 



dx° 



are the coordinate basis connection coefficients, and the 
C transformations are given by Eq. i|14fl and subsequent 
text, Eqs. (|16I2U[) and subsequent text, and Eq. IL'l'l) and 
subsequent text. We note the following combination: 



ul, v \ 

= ( ri A^V) afcc + (r^pV) Cl02C3 (37) 

The first term involves derivatives of the metric functions 
a, b, and c: 



(rW/) a . 



be 



A 3 + B 3 + C 3 , 



where 



A 3 = ( A-'j/r - A ' | //" 



p 2 da 
a dx 1 
p 3 db 



^=(aV-aV)^0 

C 3 = (\- ! ,jr - A',/- 



■] p3 dc 

c dx 2 



(38) 

(39) 
(40) 
(41) 



See Tabled for explicit values of the metric derivatives in 
various coordinate systems. The second term of Eq. I|37|) 



involves derivatives of the velocity components: 



where 



V 3 



— — -| 

dx^ 7 + 1 dx p 



= A 3 -- 

= v> 



dxP 



(42) 



7 



dvj 



Vk - 



_ dv k \ „ ' 



eW- (43) 



7 + 1 y ~ n dx^ 

We note that this contribution to Eq. I|37|) will obtain 
even in general relativity (with suitable expressions for 
the tetrad e%). 

Particularly for use in approaches employing angular 
moments, it may be worth noting that the 1 = e elements 
of the multiplication by du 1 /dp 3 simplify noticeably: 



du e 



(36) ™ A 3 



dp 3 
du 
d^ 3 
du 
dp' 3 
du e 



- - - J 2 da 

p 3 db 
~bdx [ ' 



- B 3 = 7 e (n^Vi — n^vs) 



- C 3 = je(n 3 v 2 - n 2 v 3 ) 



p 3 dc 



V 3 



enj 7 
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7+I dx^ 



(44) 
(45) 
(46) 

e%p*. (47) 



Beyond this, more explicit expressions of the 1 = i?, tp 
elements are not particularly illuminating; and even in 
computer codes it may be worthwhile, for a reason men- 
tioned in the appendix, to use the somewhat more general 
expressions in the previous paragraph. 



IV. SUMMATION 

To summarize, the nonconservative (Eq. JlJ), number- 
conservative (Eq. ©), and energy-conservative (p = 
component of Eq. ©) special relativistic transport equa- 
tions for the massless particle scalar distribution func- 
tion /, a function of global spacetime coordinates (x p ) = 

(x 1 , x 2 , x 3 , t) T and comoving- frame energy and angles 

(u l ) = (e,#, p) T — with metric functions a, b, and c suf- 
ficiently general to allow for Cartesian, spherical, and 
cylindrical coordinates — are 



(48) 



6 



_1 d_ 

abc dx^ 



1 d_ 

sin d du 1 



du 1 - 

esintf — -r (A J + 5 J + C J + y J ) / 

op? v ' 



= C[/], 



(49) 



j a_ 



_1 d_ 

mi d du 1 



du 



e sin? 



i? — r (A 3 + B 3 + C 3 + V 3 ) p°f 



(50) 



where A 3 , B 3 \ C J \ and V 3 are given by Eqs. ^39l4ip . 
(|43|1 . We remind that in all these expressions, care must 
be taken to distinguish between the momentum compo- 
nents in different frames, given by Eqs. (|f 0|) . (|23I25|) . and 
(|28l33j) . The expression dv? /dp 3 is given by Eq. (fMjl . 
Finally, 7 (given by Eq. (|20Jl ) is the Lorentz factor of 
the boost (Eqs. (|16I19|0 between the orthonormal "lab" 
and "comoving" frames. The expressions in Eq. (|44I47() 
will be convenient for use in angular moment formalisms 
(such as variable Eddington factor approaches). 

Our final equations are valid in flat spacetime. In 
Eq. H50fl . we note that r° M „ in Eq. © has vanished 
in flat spacetime, but (in general) would not vanish in 
curved spacetime. Also, p° would have to be replaced by 
the more general expression e°p l p fi . But even in curved 
spacetime, the V° term will be the same, except for a 
different tetrad e M p in Eq. I|4H|I. Finally, the major dif- 
ference in curved spacetime will be the replacement of 
Eq. (|38[) with terms arising from the more general met- 
ric. 

In the appendix we examine in detail the relationship 
between these number and lab frame energy transport 
equations, highlighting expressions that should cancel in 
a finite-differenced representation. Our finite-differenced 
representation of the special relativistic neutrino number 
transport equation and the results of core-collapse simu- 
lations using it will be reported separately, but for now 
we can mention some salient points regarding our ap- 
proach to finite differencing. First, because we must be 
able to accurately reproduce neutrino diffusion inside the 
nascent neutron star while retaining stable free stream- 
ing outside, our spatial differencing interpolates between 
centered differencing (good for diffusion) and "upwind" 
differencing (stable free streaming), according to the neu- 
trino mean free path 0, 0] . Second, also important for 
accuracy in the diffusion limit, certain terms cancel be- 
tween the space and momentum space divergences in the 
homogeneous and isotropic limit; differencings of A J , B J , 
and C 3 are chosen that respect this 0,0. Third, the 
conservative formulation — with its divergences in space 
and momentum space — invites a conservative differenc- 
ing, patterned after the mathematical definition of the 
divergence: a sum over cell faces, divided by the cell 
volume. Global conservation then follows naturally, fol- 
lowing the discrete version of the divergence theorem. A 
discrete volume integral is obtained by multiplying the 
divergence in each cell by the cell volume; the contribu- 
tions on faces shared by adjacent cells then cancel ex- 
actly. Fourth, because neutrino transport in supernovae 



requires both number and energy conservation, we dis- 
cretize the number-conservative Eq. I|49|l . and work out 
algebraically the differencings required to ensure consis- 
tency with lab-frame energy conservation as encoded in 
Eq. (|50|l . This involves, for example, the differencing 
of velocity gradients in V 3 beingtied to the differencing 
chosen for the space divergence 

While our parochial concern is neutrino transport in 
supernovae, relativistic photon transport is important in 
many astrophysical systems. Compton scattering is an- 
alyzed in detail in Refs. H3,|2^|, where it is concluded 
that 0{v/c) transfer equations are inadequate even for 
electron bulk velocities v ~ 0.1c. This work was ap- 
plied to spherically symmetric accretion |28| . Obvious 
examples in which multiple space dimensions are impor- 
tant are gamma-ray bursts and "superluminal" galactic 
jets. A ray-tracing approach (as opposed to the grid- 
based discretization we intend to employ) in a back- 
ground stationary Kerr metric showed that both kine- 
matic and spacetime curvature effects are important in 
photon transfer in accretion disks and tori surrounding 
black holes [2j|. Another example is photon transport 
in supernovae, necessary for the interpretation of spec- 
tra and lightcurves. It has been shown in spherically 
symmetric simulations that all special relativistic terms 
must be included to accurately model the spectra of the 
modestly relativistic (v ~ 0.1c) ejecta of supernovae 
|30|. This fully relativistic formulation has been rou- 
tinely used to evaluate both spherically symmetric explo- 
sion models and observed objects of all supernova types 
[3TI I32I I33I l34| . Observational indications that excep- 
tional supernovae may produce ejecta that is both rela- 
tivistic and jet-like |35j suggest that similarly relativistic 
completeness in multidimensional modeling of these ob- 
jects' spectra and lightcurves will be necessary. 



APPENDIX 

It is apparent from the right-hand sides of Eqs. 149|) 
and (|50fl that the latter is simply p° times the former, but 
it is not obvious how this works out on the left-hand side. 
Because p° is inside the spacetime and momentum space 
divergences in Eq. (|50|l . the "extra" terms generated by 
pulling p° inside the derivatives of Eq. (|49|l must cancel. 
Here we explore this cancellation in some detail. 

Consider first the "extra" term arising from the space- 
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time divergence: 



p 5 d 
abc dx^ 

where 
Es 



(abce^ppPf) 



_L d_ 

abc dx^ 



(A.l) 



- dp° 



dx^ 

7^1^37 + + W 



(A.2) 
,(A.3) 



and Eq. (|29|) has been used in the last line. 

We begin consideration of the "extra" term resulting 
from the momentum divergence by noting that care must 
be taken in the calculation of momentum derivatives of 
p . For this purpose, uncritical use of the expression in 
Eq. (|29[) is misleading. Recall that p° is the result of a 



Lorentz boost from the comoving frame; and that while 
pV- is a 4-vector, because it is on the mass shell (of a 
massless particle in the present case), p° — e is to be 
considered a function of the p % : 



P° = 



(A.4) 



= AV + A%^l (pi)' + ( P 2) 2 + ^i) 2 . (A.5) 



Hence 



dp = d^<¥ = 1 / eA o. + A o. \ M 
du l dp 1 du l e \ 1 V du l 



(A.6) 



is the needed derivative of this component. 

Now we may contemplate the "extra" term resulting 
from the momentum divergence: 



p° d 
e sin du l 

where 



esin-i? — (A J + B 3 + C 3 + V J ) f 
dpi v ' ' 



___L d_ 

e sin •& du l 



esintf — (A 3 + B 3 + C 3 + V J ) p / 



£m, (A.7) 



Em = ~i$( a3+bS+cS+v3 )^ 



-i (e A 5 ,- + a ',,/,,) + /;•' + r' + v 5 ) /. 



(A.8) 
(A.9) 



In obtaining the last line, we used Eq. (|A.6|) . which re- 
sulted in a delta function as a consequence of the multi- 
plication of the momentum space Jacobian du l /dp 3 with 
its inverse dp l /du l . In the last sentence of Sec. IIIII 
we referred to the possible utility of retaining even in 
computer codes — our expressions involving the matrix 
du 1 /dp 3 , and we are now in a position to see why: The 
finite-difference representation of dp 1 /du l the right-hand 
side of Eq. (|A.6(1 will be determined by a choice of the 
finite-difference representation of the momentum space 
divergence; hence the finite-difference representation of 
du 1 /dp 3 can be chosen to make it precisely the matrix 
inverse of the (foreordained) finite-difference representa- 
tion of dp 1 /du % . 

We expect the cancellation E$ + Em = in accordance 
with Eqs. (|49|l and l]5Up. and because Eg contains no 
derivatives of the metric functions a, b, c, there must be 
"internal cancellations" in Em such that 

(e A°j + A° 6 pj) (A 3 + B 3 + C 3 ) = 0, (A. 10) 

as we now verify. From Eqs. (|39I41I) we see that each 
of the terms arising from A 3 , B 3 , and C 3 has an overall 



multiplicative factor of the form 

(e A°j + A° 6 Pi ) (.V. // - A* 5 //) . ( A. 1 1 ) 

We note that (for example) 

A 5 3 A\ = A^A^-A^A 6 , (A.12) 
= -A° 6 A° f . (A.13) 

(The first term in Eq. (|A.12|I vanishes, being equal to 
<5V) Hence 

'•V', (AW - As ,/) = -e k\ (A ',// - A',,.,/) . 

(A.14) 

Furthermore, we note that (for example) 

PjA\ = p A A^-p 5 A 5 , (A. 15) 

= p* + eA 6 f (A.16) 

(recall that p^ = rj^^pi 1 = — p°), so that 
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"„/>, (AW - As //) = A",, ,>,,,>') + < (A 8 ,/ - A",,.//) 

= £ a 5 6 (aV-a 6 ^ 



(A.17) 
(A.18) 

(The first term in Eq. I|A.17|) vanishes because pi — p l in an orthonormal basis.) Taken together, Eqs. I|A.14J| and 
(|A.18|I imply that Eq. (|A.11J| vanishes, proving Eq. I|A.10|) as desired. 
The "extra" term Em is now reduced to 



Em — 



A 5 j + A 5 5 p/) V'f 



2 dvj jVj 



dx^ 7+1 dx^ 7+1 



dvj 



Vh Vn 

dx» 3 dx^ 



>ik 



(A.19) 
(A.20) 



where we have used Eqs. (|18fl . I|19|) . and (|23I25(I to rewrite the leading factor, and Eq. (|43|l to substitute for V' J . We 
denote the quantity in square brackets by [...]. The last two terms of [. . .] don't contribute to the combination hj [...], 
because the contraction of the symmetric combination hjhk with an expression antisymmetric in j and k vanishes; 
all that remains is 



7 rij 



dvj_ rhjVj 07 
g x n 7 + 1 dx^ 



Using the identities VjVj = (7 — 1)(7 + l)/7 2 and Vj(dvj/dx^) — 7 3 (d^/dx fi ), we find that 

1,7-1, fi k v k \ 97 _ dv k 

7 7 7 + 1 / aa;' 1 ox'' 

\ ^7 _ ( 2 _ } 9^ 



Plugging the sum of Eqs. (|A.21jl and (|A.23jl into Eq. (|A.20|) yields 

E M = -ee»-^f 



7^ + (!+»*«*) fe? 



(A.21) 

(A.22) 
(A.23) 

(A.24) 



Comparison with Eq. (|A.3|I for i?s shows that the desired cancellation E$ + Em = is verified. 
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